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ABSTRACT 

We report the detection of CO(J=l— »0) emission in the strongly lensed high-redshift 
quasars IRAS F10214+4724 0=2.286), the Cloverleaf 0=2.558), RXJ0911+0551 0=2.796), 
SMM J04135+10277 0=2.846), and MG 0751+2716 0=3-200), using the Expanded Very Large Ar- 
ray and the Green Bank Telescope. We report lensing-corrected CO(J=l— >0) line luminosities of 
L' co = 0.34-18.4xl0 10 K kms _1 pc 2 and total molecular gas masses of M(H 2 ) = 0.27— 14. 7x 10 10 Mq 
for the sources in our sample. Based on CO line ratios relative to previously reported observations in 
J>3 rotational transitions and line excitation modeling, we find that the CO(J=1^0) line strengths 
in our targets are consistent with single, highly-excited gas components with constant brightness tem- 
perature up to mid- J levels. We thus do not find any evidence for luminous extended, low excitation, 
low surface brightness molecular gas components. These properties are comparable to those found 
in z>4 quasars with existing CO(J=l— >0) observations. These findings stand in contrast to recent 
CO(J=l— >0) observations of z~2-4 submillimeter galaxies (SMGs), which have lower CO excitation 
and show evidence for multiple excitation components, including some low-excitation gas. These find- 
ings are consistent with the picture that gas-rich quasars and SMGs represent different stages in the 
early evolution of massive galaxies. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Investigations of the molecular and dusty interstellar 
medium (ISM) in high redshift galaxies are of key im- 
portance to studies of galaxy evolution at early cosmic 
times, as the ISM provides the material that fuels star 
formation and stellar mass assembly. It is particularly 
interesting to better understand the ISM properties of 
galaxies that host both intense star formation and a lu- 
minous active galactic nucleus (AGN), as this enables 
simultaneous investigations of supermassive black hole 
and stellar bulge growth in galaxies. A key cosmic epoch 
for such studies is the redshift range 2<z<3, where both 
cosmic star formation and AGN activity peak, and thus, 
where most of the growth of stellar an d bla ck hole mass 
in galaxies occurs (e.g., Magnelli et al. 120091: Richards et 
al 



gaiaxi 
2006) 



1 Astronomy Department, California Institute of Technology, 
MC 249-17, 1200 East California Boulevard, Pasadena, CA 
91125, USA; dr@caltech.edu 

2 National Radio Astronomy Observatory, PO Box O, Socorro, 
NM 87801, USA 

3 National Radio Astronomy Observatory, P.O. Box 2. Green 
Bank, WV 24944, USA 

4 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D- 
69117 Heidelberg, Germany 

5 Department of Astronomy, University of Maryland, College 
Park, MD 20742-2421, USA 

6 Department of Physics and Astronomy, Rutgers, the State 
University of New Jersey, 136 Frclinghuysen Road, Piscataway, 
NJ 08854-8019, USA 

7 European Southern Observatory, Alonso de Cordova 3107, 
Vitacura, Casilla 19001, Santiago 19, Chile 

8 National Radio Astronomy Observatory, 520 Edgemont 
Road, Charlottesville, VA 22903-2475, USA 

9 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigcl 
69, D-53121 Bonn, Germany 



Some of the most remarkable distant AGN-starburst 
galaxies are gas-rich, far-infrared luminous quasars at 
z>2. To date, the molecular ISM of 34 high-z quasars 
has been detected through emission from rotational lines 
of CO, 14 of which are gravitationally lensed (s ee rev iew 
by Solomon & Vanden Bout|2005j, and Riechers HoTJ for 
a recent summary). Given the sensitivity of past obser- 
vatories, gravitational lensing is key to probe down to 
intrinsically fainter (and thus, more common) systems. 
One key aspect of these studies is that most detections 
were obtained in mid-J (i.e., J>3) CO transitions. Only 
three of the 34 quasars were observed in the ground- 
state CO(J= 1->0) line (e.g., Carilli et al. 120021: Riechers 
et al. I2006a|) . It thus remains a possibility that many 
of these studies are biased toward the highly excited gas 
that does not necessarily trace the entire molecular gas 
reservoir seen in CO( J—l— >0). 

To overcome the limitations of previous studies, we 
have initiated a systematic study of the CO(J=l— >0) 
content of high-z quasars and other galaxy populations 
wit h the Expanded Very Large Array (EVLA; Perley et 
al. 1201 If) and the 100 m Robert C. Byrd Green Bank 
Telescope (GBT). In this Letter, we report the detec- 
tion of CO( J=l— >0) emission toward five strongly lensed 
z>2 quasars, using the EVLA and the GBT. We use 
a concordance, flat ACDM cosmology throughout, with 
H n = 71 kms-iM pc" 1 , ft M =0.27, and S1 A =0.73 (Spergel 
etal. l200ai2007h . 

2. OBSERVATIONS 
2.1. EVLA 

We observed the CO(J=1^0) (z/ rest = 115.2712GHz) 
emission line toward IRAS F10214+4724 0=2.286), the 



Riechers et al. 



0.0 


[mjy beam] 
0.5 1.0 1.5 


1 1 . T 


~CO(J= 


1-0) in p1 021 4+4724 (z=2.286) 












2 .7"x2.-3 n :(D Configuration 



47 09 25 
f 20 

CM 
-> 

Z 15 

g 
t- 

< 10 

o 

S 05 
00 
08 55 

10 24 36.0 35.5 35.0 34.5 34.0 33.5 33.0 
RIGHT ASCENSION (J2000) 




^^^^^^^^^^ 


CO(J 


I I I < ~! 

=1-0) in RX J091 1 +055fr (z=2.796) 




i 






m 




; « 






-0 




O 




3.3"x2.4" (D configuration) 



27 16 55 
50 
45 
40 
35 
30 
25 
20 
15 
10 



5 


10 


I i ~Tj 


~CO(J=1-0) in MG 0751+2716 (z=3.200)~ 


_ +continuum 

w 






co(jji-o) m 








4 <f I 


X)| 3.1"x2.4"(Dconfig) 





1415 47.0 46.6 46.2 45.8 45.. 



09 11 28.5 28.0 



07 51 43.0 42.0 



Fig. 1.— EVLA CO(J=l-s>0) maps of IRAS F10214+4724 (2=2.286), the Cloverleaf (2=2.558), RX J0911+0551 (z=2.796), and 
MG 0751+2716 (z=3.200). The resolution is indicated at the bottom of each panel. For MG 0751+2716, the inset shows the CO line 
emission after subtracting the underlying continuum emission. Contours are in steps of 1<t=125, 115, 75, and 330 /^Jy beam" 1 for each 
source, respectively, starting at ±2cr (except the CO+continuum map of MG 0751+2716, where contours are in steps of 4<r). The maps are 
averaged over 160, 648, 185, and 306 kms -1 (18.75, 70, 18.75, and 28 MHz), respectively. 



Cloverleaf (z=2.558), RXJ0911+0551 (^=2.796), and 
MG 0751+2716 0=3.200), using the EVLA. At these 
rcdshifts, all lines are shifted to the Ka band (0.9 cm; 
see Table Q] for redshifted line frequencies). Observations 
were carried out under good weather conditions in six 
D array tracks between 2009 October 26 and Decem- 
ber 09, and on 2010 July 17 and 18. This resulted in 
5.0, 1.0, 12.0, and 1.0 hr (2.1, 0.6, 5.1, and 0.6 hr) total 
(on-source) observing time for IRAS F10214+4724, the 
Cloverleaf, RX J0911+0551, and MG 0751+2716, respec- 
tively. For IRASF10214+4724 and RX J0911+0551, an 
additional 0.9 and 2.3 hr on source were spent on separate 
continuum settings. The nearby quasars J0958+4725, 
J1415+1320, J0909+0121, and J0748+2400 were ob- 
served every 3.5 to 7 minutes for pointing, secondary am- 
plitude and phase calibration. For primary flux calibra- 
tion, the standard calibrators 3C286 and 3C147 were ob- 
served, leading to a calibration that is accurate within 
-10%. 

Observations for the Cloverleaf and MG 0751+2716 
were carried out with the WIDAR correlator, using two 
intermediate frequencies (IFs) of 128 MHz (dual polar- 
ization) each at 2 MHz resolution. For the Cloverleaf, 
the two IFs were overlapped by two channels, centered 
on the CO line, yielding 252 MHz contiguous bandwidth. 
For MG 0751+2716, one IF was centered on the CO line, 
and the second IF was centered on the continuum at 
32.046 GHz. 

Observations for IRAS F10214+4724 and 
RX J0911+0551 (which have narrow CO lines) were 
carried out with the previous generation correlator, with 
two 21.875 MHz (dual polarization) IFs at 3.125 MHz 
resolution. For IRAS F10214+4724, both IFs were 
centered on the CO line, yielding 43.75 MHz contiguous 
bandwidth. For RX J0911+0551, one IF was centered 
on the CO line, and the second IF was centered on 
the continuum at 34.3173 GHz. For these two sources, 
one third of the on-source time was spent to observe a 
second frequency setting with two 50 MHz continuum 
IFs offset by ±150 MHz from the CO lines, yielding 
more sensitive constraints on the continuum emission. 

For data reduction and analysis, the AIPS pack- 
age was used. All data were mapped using 'natural' 
weighting. The data result in final rms noise levels of 
125, 115, 75, and 330 /iJy beam" 1 over 160, 648, 185, 
and 306 kins" 1 (18.75, 70, 18.75, and 28 MHz) widths 
for IRASF10214+4724, the Cloverleaf, RX J0911+0551, 



and MG 0751+2716, respectively. Maps of the velocity- 
integrated CO J=l— >-0 line emission yield synthesized 
clean beam sizes of 2.7"x2.3", 3.7"x2.1". 3.3"x2.4", and 
3.1"x2.4". 

2.2. GBT 

We observed the CO(J=l— >0) emission line toward 
IRASF10214+4724, the Cloverleaf, SMM J04135+10277 
(z=2.846), and MG 0751+2716, using the facility Ka 
band receiver on the GBT with the digital Spectrom- 
eter backend (GBT/S) and the Zpectrometer analog lag 
cross-correlation spectrometer backend (GBT/Z; instan- 
taneously covering 25.6-36.1 GHz; Harris et al. 120071 : 
see Table [1]). This yields a typical beam size of —23". 
GBT/S observations were carried out under acceptable 
to good weather conditions during 14 sessions between 
2007 October 10 and 2008 October 12, yielding typi- 
cal system temperatures of T sys =42-57 K (with higher 
T sys toward higher observing frequencies due to the at- 
mosphere). This resulted in 4.3, 3.7, 2.5, and 8.7 hr 
on-source observing time for IRAS F10214+4724, the 
Cloverleaf, SMM J04135+10277, and MG 0751+2716, re- 
spectively. In addition, the Cloverleaf was observed with 
the GBT/Z for another 6 sessions between 2008 March 02 
and 2009 March 12. This resulted in 14.7hr of observing 
time, about half of which was spent on the Cloverleaf. 
The GBT/S was configured with 800 MHz bandwidth, 
yielding a spectral resolution of 391kHz. The GBT/Z 
samples its 10.5 GHz bandwidth with 8 MHz spectral 
channels. Its instrumental spectral response is nearly 
a sine function with an FWHM of 20 MHz, i.e., individ- 
ual 8 MHz spectral channels are not statistically indepen- 
dent. However, the line width correction for the instru- 
mental response for spectral lines with intrinsic Gaussian 
FWHM of >30MHz (-300 kms^ 1 ) is minor. 

Subreflector beam switching was used every 10 s for 
all observations (GBT/S and GBT/Z) to observe the 
sources alternately with the receiver's two beams, with 
the off-source beam monitoring the sky background 
in parallel. To remove continuum fluxes and atmo- 
spheric/instrumental effects, low-order polynomials were 
fitted to the the spectral baselines of the calibrated, time- 
averaged GBT/S observations. For the GBT/Z observa- 
tions, a nearby second, faint source was observed with the 
same subreflector switching pattern, alternating between 
targets with 8 min cycles. Residual structure from optical 
beam imbalance in the "source-sky" difference spectra of 
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TABLE 1 

Observed CO(J=l— »-0) line parameters. 



Source 


z CO 


Ml 


^obs 

[GHz] 


S v 
\vc\ Tvl 


AVpwHM 

[ km s ] 


I CO 
f Tv km s 1 


j 1 a 

[10 9 K kms -1 pc 2 ] 


Telescope 


F10214+4724 


2.2853 ± 0.0001 
2.2856 ± 0.0001 

2.2854 ± 0.0001 


17 


35.0795 


2.42 ± 0.37 
1.73 ± 0.22 


169 ± 39 
184 ± 29 


0.434 ± 0.047 
0.337 ± 0.045 
0.383 ± 0.032 


5.77±0.49 


EVLA 
GBT/S 
Combined 


Cloverleaf 


2.5564 ± 0.0004 
2.5578 ± 0.0001 

2.5574 ± 0.0001 

2.5575 ± 0.0001 


11 


32.3992 


2.78 ± 0.40 
2.72 ± 0.12 
3.14 ± 0.11 


468 ± 94 
470 ± 27 
422 ± 19 


1.378 ± 0.250 
1.358 ± 0.065 
1.406 ± 0.053 
1.387 ± 0.040 


39.3±1.1 


EVLA 

GBT/S 
GBT/Z 
Combined 


J0911+0551 


2.7961 ± 0.0001 


22 


30.3665 


1.75 ± 0.22 


111 ± 19 


0.205 ± 0.029 


3.39±0.48 


EVLA 


J04135+10277 


2.8470 ± 0.0004 


1.3 


29.9717 


1.20 ± 0.15 


505 ± 75 


0.644 ± 0.082 


184±23 


GBT/S 


MG 0751+2716 


3.1984 ± 0.0006 
3.1995 ± 0.0004 
3.1990 ± 0.0003 


16 


27.4455 


1.61 ± 0.34 
1.47 ± 0.24 


290 ± 62 
354 ± 72 


0.494 ± 0.105 
0.550 ± 0.095 
0.525 ± 0.070 


14.9±2.0 


EVLA 

GBT/S 
Combined 



Note. — fi^: lensing magnification factor (see Riechers |2011| , and references therein). GBT/S: Measured with the GBT digital Spectrometer 
backend. GBT/Z: Measured with the GBT wideband Zpcctromctcr backend. 



a Corrected for magnification due to gravitational lensing. 
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Fig. 2.— Overlay of CO(J=1->0) and CO(J=7->-6) emission (left; CO J=7->6 map from Alloin et al. 1997 ). and independent EVLA, 
GBT/Spectrometer, and GBT/Zpectrometer CO(J=l— >0) spectra (middle left to right) in the Cloverleaf (z=2.558). The spectra (his- 
tograms) are shown at 4, 3.906, and 8 MHz resolution (37, 36, and 74 kms -1 ). The solid curves indicate Gaussian fits to the spectra. The 
GBT/Z channels are not statistically independent on scales below 20 MHz (2.5 channels). 

the two targets was then removed by differencing the re- 
sulting spectra of both sources (the second source was not 
detected) . This strategy yields a flat baseline (offset from 
zero flux by the difference of the two source's continua) 
without standard polynomial baseline removal. For all 
observations, several nearby quasars were targeted reg- 
ularly to monitor telescope pointing and gain stability. 
Passband gains and absolute fluxes for all observations 
were determined from spectra of 3C286 and 3C147, yield- 
ing 10%— 15% calibration accuracy. 
All da ta were processed with GBTIDL (Marganian et 



al.[2b06), using standard recipes and/or the Z pectro m- 
eter's data reduction pipeline (see Harris et al. 120101 for 
GBT/Z details). After processing, the GBT/S data were 
re-binned to 3.906 MHz (33-43 kms -1 ) for further anal- 
ysis. 

3. RESULTS 

We have detected strong CO( J—l— >0) emission toward 
IRASF10214+4724, the Cloverleaf, RX J0911+0551, 
SMMJ04135+10277, and MG 0751+2716. We spatially 
resolve the CO( J=l->0) emission in IRAS F10214+4724 
and the Cloverleaf, but do not resolve the CO(J=l— ?-0) 
emission in RX J0911+0551 and MG 0751+2716 (Fig. [I]). 
The measured sizes and upper limits are consistent with 
the lens configurations of these targets if the CO reser- 
voirs have intrinsic sizes of a few kpc. Accounting for 
beam convolution, the gas distribution in the Clover- 
leaf is consistent with that observed in higher- J lines 



(e.g., Alloin et al. 1997), but may be somewhat more ex- 
tended (Fig. [21 left). We detect rest-frame 2.6mm contin- 
uum emission toward the Cloverleaf and MG 0751+2716 
at 386+107 /iJy and 17.4+0.7 mJy strength. We also 
detect rest-frame 2.2 mm continuum emission toward 
MG 0751+2716 at 13.9+0.6 mJy strength. We do not 
detect rest-frame 2.6 mm continuum emission toward 
IRAS F10214+4724 and RX J0911+0551 down to 3ct lim- 
its of 270 and 170 /iJy. The continuum emission in 
MG 0751+2716 is consistent with non-ther mal e mission 
from its radio-l oud AGN (e.g., Lehar et al. 1 19971 Riech- 
ers et al. |2006b[) . The continuum detection and limits for 
the other targets are consistent with thermal and non- 
thermal emission associated wi th sta r formati on in their 
host galaxi es (e.g., Weifi et al. 120031 : Ao et al. 120081 : Wu 
et al. 12009(1 . 

The CO(J=l— >0) line profiles of all sources were fit- 
ted with single-component Gaussians, which yields good 
fits to all data (Figs. [2] and [3]). The three independent 
measurements of the CO( J=l— >0) line in the Cloverleaf 
(EVLA, GBT/S, and GBT/Z) yield consistent line pro- 
files and intensities within the relative uncertainties, sug- 
gesting that the interferometric observations recover the 
full flux, and that the calibration of the single-dish ob- 
servations is reliable to high accuracy. Due to the lim- 
ited bandwidth of the observations of RX J0911+0551, 
the line is only marginally covered in velocity. How- 
ever, the fitted line FWHM of 111+19 kms -1 is in ex- 
cellent agreement with the ~110 kms -1 measured in the 
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Fig. 3.— EVLA and GBT/Spectrometer spectra of CO(J=l-s>0) emission toward IRAS F10214+4724 (2=2.286), RX J091 1+0551 
(2=2.796), SMMJ04135+ 10277 (2=2.846), and MG0751+2716 (2=3.200). The spectra (histograms) are shown at 3.906MHz resolu- 
tion (33-43 kms -1 ), except RX J0911+0551, which is shown at 3.125MHz (31 kms" 1 ) resolution. The solid curves indicate Gaussian fits 
to the spectra. The blanked-out region in the spectrum of MG 0751+2716 masks a resonance from the GBT Ka band receiver's feeds. 



CO(J=3^2) to CO(J=9^8) lines (A. Weifi et al., in 
prep.)E3 The GBT/S measurement of the CO(J=1^0) 
line in MG 0751+2716 shows an instrumental resonance 
feature —800 kms -1 bluewards from the line center, 
which is excluded from the spectral baseline fit, and does 
not appear to affect the line flux measurement. The 
CO( J=l—>0) line peak flux densities, widths, integrated 
intensities, and centroid redshifts for all sources and mea- 
surements are detailed in Table [TJ Two or three inde- 
pendent CO(J=l— ?>0) measurements were obtained for 
IRASF10214+4724, the Cloverleaf, and MG 0751+2716. 
The median CO redshifts and line intensities for these 
sources are given in Table [TJ and are adopted for all fur- 
ther analysis. All five targets are gravitationally lensed. 
Based on the CO(J=l— >-0) line intensities and lensing 
magn ification factors from the literature (see Riechers 
120111 and references therein), we derive lensing-corrected 



CO(J=l— >-0) line luminosities of 

18.4 x 10 10 K kms' 

Our CO(J=1^0) 
CO(J=3^2)/CO(J=1^0) line 



^pc 2 (Tab. [JJ). 



ji 

^CO(l-O) " 

fluxes 
brightness 



CO(3-2) /T ,CO(1-0)_ 



0.34 

imply 
temper- 

1.00±0.10, 



ature ratio£3 of r 31 =T^~*> /T h 
1.06±0.03, -0.95, 0.93±0.25, and 0.97±0.17 for 
IRASF10214+4724, the Cloverleaf, RX J0911+0551, 
SMMJ04135+10277, and MG 0751+2716, respective! 
(CO J=3 — >2 m easurements are from Ao et al. 1200 
Weifi et al. [20031 A. Weifi et al., in prep., Hainline et al. 
12004 and Alloin et al. 120071 respectively). These values 
are fully consistent with thermalized gas excitation (i.e., 
r3i=l) in all targets. 

4. ANALYSIS 

4.1. Line Excitation Modeling 

The two best-studied targets in our sample are 
LRAS F10214+4724 a nd th e Cloverleaf (e.g. Ao et al. 
[2001 Barvainis et al. [19971 Weifi et al. 120031: Bradford 
et al. [2009I ). Based on CO(J=3^2) to CO(J=9^8) 
observations in the literature and our CO(J=l— >Q) de- 
tections, we can constrain the line radiative transfer 
through Large Velocity Gradient (LVG) models, treat- 
ing the gas kinetic temperature and density as free pa- 
rameters. For all calculations, the H2 ortho-to-para ra- 

10 These measurements revise the C O(J=3— >2) line flux and 
width measured by Hainline et al. 20041 down by ~40% and by a 
factor of ~3, respectively. 

11 The error bars for r^i are derived from the statistical uncer- 
tainties of the CO(J=l— ¥0) and CO(J=3— >2) measurements. 



tio was fixed to 3:1, the cosmic microwave background 
temperature was fixed to 8.95 and 9.69 K (at 2=2.286 
and 2.558), and the Flower (|200lD CO collision rates 
were used. For consiste ncy w ith previous mod eling o f 
both sources (Ao et al. 120081: Riechers et al. 1201 laf ). 
we adopted a CO abundance per velocity gradient of 
\CO]/(dv/d r) = 1 x lQ- 5 pc(km s- 1 )~ 1 (e.g., Weifi etal. 
[20051 [2001 Riechers et al. l2Q06ah . 

Observations of both targets are fit very well by single- 
component models (except the CO J=5—>4 flux in the 
Cloverleaf, which is likely too low due to calibration is- 
sues related to the restricte d ba ndwidth of these early 
observations; Barvainis et al. l!997t see also A. Weifi et al., 
in prep.). For IRAS F10214+4724, we fit a representative 
model with a kinetic gas temperature of Tki n =60K and 
a gas density of /? gas =10 3 8 cm" 3 , yielding a moderate 
optical depth of rco(i-o)— 1-6 (Fig. [4ja and c). For the 
Cloverleaf, we fit a representative model with Tki n =50K 
and /5 g as=10 4 ' 5 cm" 3 , yielding a high optical depth of 
t CO (i-o)=8.9 (Fig.[4p and d). While not formally exclud- 
ing the presence of some colder gas (given the remaining 
uncertainties), these models are consistent with previous 
fits b ased on the C O J>3 transitions al one (A o et al. 
12001 Bradford et al. l20"09l Riechers et al. l2011afl . They 
are also consistent with what is found for z>A quas ars ob- 
served in C O( J= l— >0) emission (Riechers et al. I2006at 
Weifi et al. |2007[ ). This finding is in agreement with 
the CO( J=l— >Q) emission in high-z quasars being associ- 
ated with optically thick emission from the highly excited 
molecular gas detected in high-J CO transitions, with- 
out any evidence for significant additional low-excitation 
gas components. 

4.2. Molecular Gas Masses 

The CO(J=l— >Q) emission in our targets appears to 
be associated with highly excited gas that has physical 
properties consistent with those found in the nuclei of 
nearby ultra- luminous infrared galaxies (ULIRGs). We 
thus calculate their gas masses assuming a ULIRG-likc 
conversion factor of aco=0.8 M Q (K kms -1 pc 2 ) -1 from 
^CO(i-o) to ^gas (Downes & Solomon ll998l ). This yields 
lensing-corrected molecular gas masses of M gas =0.46, 
3.1, 0.27, 14.7, andl.2xlO lo Af for IRAS F10214+4724, 
the Cloverleaf, RX J0911+0551, SMM J04135+10277, 
and MG 0751+2716, respectively. SMM J04135+10277 is 
now revealed as the high-z quasar with the highest M gas 
currently known. The gas masses of all systems lie within 
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Fig. 4. — CO excitation ladders (points) and LVG models (lines) for IRAS F10214+4724 (panels a and c) and the Cloverleaf (panels b 
and d), scaled by CO intensity (left two panels) and normalized brightness temperature Tb (right two panels). The insets show zoo med-in 
versions of the model fits close to the CO(J=1->0) line. The C O dat a for IRAS F10214+4724 ( circle : this work; squares: Ao et al. 120081) 
and the Cloverleaf (circle: this work; triangles: Barvainis et al. 119971 : filled square: Weifi et al. 20031; filled circles: Bradford et al. 2009; 
error bars of literature data have been conservatively extended to include 10%-15% calibration uncertainties) are fit well by single, highly 
excited gas c ompon ents with kine tic tem peratures of Ti c : n =60 and 50K and gas densities of p ga s=10 3 ' 8 and 10 4 ' 5 cm -3 , respectively (see 
also Ao et al. [200g ; Riechers et al. l2011aT) . 

the range found for other high-z quasars (mostly deter- 
mined based on observations of J>3 CO transitions), 
with IRASF10214+4724 and RXJ0911+0551 being sit- 
uated at the low end of the observed range (which can 
currently only be investigated with the aid of strong grav- 
itational lensing; e.g., Riechers 1201 

5. DISCUSSION AND CONCLUSIONS 

We have detected strong CO( J=l— >0) emission toward 
IRASF10214+4724 (z=2.286), the Cloverleaf (z=2.558), 
RXJ0911+0551 (3=2.796), SMM J04135+10277 
(2=2.846), and MG 0751+2716 (2=3.200). Our 
EVLA observations have spatially resolved the emission 
toward IRAS F10214+4724 and the Cloverleaf (Fig. QJ. 
We find line brightness temperatures consistent with 
those measured in mid- J CO lines. Excitation modeling 
indicates that the CO(J=l— >Q) emission is associated 
with the warm, highly excited gas in the star-forming 
regions in the host galaxies that is also seen in the 
higher-J CO lines. This result suggests that mid- J 
CO lines are good indicators of the total amount of 
molecular gas in gas-rich high redshift quasar host 
galaxies, consistent with previous findings bas ed on a 
smaller sample of z>4 quasars (Riechers et al. I2006al ). 
and extending these studies to the peak epoch of cosmic 
star formation and AGN activity. In contrast, recent 
studies of CO(J = l— >Q) emission in z>2 SMGs suggest 
that many of these galaxies have substantial amounts 
of lo w-exci tation gas (e. g., Ha inline et al. 120061: Carilli 
et al. I2010t H arris et al. l2010t Ivi son et al. 120101 120111: 
Riechers et al. I2010L l2011bl l2011dh Such a difference is 
in agreement with the picture that gas-rich quasars and 
SMGs represent different stages in the early evolution 
of massive galaxies. Such an observational finding 
would be consistent with a high redshift analogue of the 
ULIR G-qua sar transition scenario proposed by Sanders 
et al. lfl98l) . 

In this picture, it may be expected to find some sources 
with gas properties that overlap with both populations, 
i.e., SMGs with molecular gas properties closer to those 



of the quasars studied here. The recently identified 
z=2.957 SMG HLSW-01 may be an example of such 
"overlap" sources. Most of its gas properties are con- 
sistent with those of "typical" SMGs (i.e., its gas mass, 
gas mass fraction, gas depletion timescale, star forma- 
tion efficiency, specific star f ormati on rate, and dynami- 
cal structure; Riechers et al. 1201 ldh . but it shows a high 
r 3 i=0.95±0.10 (i.e., consistent with 1), an d an overall 
high CO excitation (Riechers et al. 1201 ldt Scott et al. 
2011). Its high radio luminosity, dust temperature, and 
CO excitation may suggest the prese nce of a luminous, 
per haps o bscured AGN (Con ley et al. 120111 : Riechers et 
al. 1201 ldt Scott et al. 120111) . However, based on CO 
studies at high z to date, systems like HLSW-01 appear 
to be rare among SMGs. 

These results demonstrate that observations of low- 
order CO transitions in high redshift galaxies with the 
EVLA will be of enormous utility to establish the con- 
text for studies of highly excited CO emission with the 
Atacama Large (sub)Millimeter Array (ALMA). Com- 
plementary studies with the EVLA and ALMA will be 
key to distinguishing different high-z galaxy populations 
based on the physical properties of the molecular gas in 
their star-forming environments. 
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